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We study the conductance of an interconnect between two graphene leads formed by a single-atom 
carbon chain. Its dependence on the chemical potential and the number of atoms in the chain is 
qualitatively different from that in the case of normal metal leads. Electron transport proceeds 
via narrow resonant states in the wire. The latter arise due to strong reflection at the junctions 
between the chain and the leads, which is caused by the small density of states in the leads at low 
energy. The energy dependence of the transmission coefficient near resonance is asymmetric and 
acquires a universal form at small energies. We find that in the case of leads with the zigzag edges 
the dispersion of the edge states has a significant effect on the device conductance. 
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I. INTRODUCTION 

The efforts to miniaturize electronic devices have long 
motivated studies of electron transport in molecular and 
atomic scale devices. Creating reliable electrical contacts 
with the molecule presents a major challenge in molecular 
electronics. Until recently, in molecular electronic devices 
the molecule was typically att ached bet ween two normal 
metal, such as gold, electrodes ^ 1 2 1 3 1 4 1 5 1 6 1 . 

Carbon-based conductors have long been expected to 
be promising components of electronic devices-^. Apart 
from bulk graphite there are also quasi-one dimensional 
(carbon nanotubes^) and two-dimensional (graphene^) 
forms, which have remarkable mechanical and electri- 
cal properties and form strong chemical bonds with each 
other. This offers the prospect of building entire elec- 
tronic devices or circuits out of carbon-based materials. 

Single-atom carbon chains (SACCs) are natural com- 
ponents of such devices. They are expected to be ideal 
one-dimensional conductors^. They covalently bond to 
other carbon materials. Formation of SACCswas conjec- 
tured to occur between carbon nanotube d 11 ! 12 ! 13 ! and in 
gaps between two graphene leads^. Formation of SACCs 
between graphene electrodes fabricated by stretching a 
graphene strip has been observed in Ref. [TS] SACC inter- 
connects between graphene leads could form a basic unit 
for integration into more complicated circuits in the fu- 
ture. Electron transport through SACCs with graphene 
contacts has not been studied theoretically. We address 
this issue in the present paper. 

Electron transpor t in SACC s with metal leads has been 
studied numeric all y! 1 * 1 6 * X 211S1 ^he conductance of an 
SACC strongly coupled to metal electrodes is of order 
of the conductance quantum and exhibits even-odd oscil- 
lations with the number of atoms in the chain, with the 
contrast ratio of order unity. 

In this paper we show that electron transport through 
an SACC interconnect between graphene leads in the 



non-interacting electron approximation is described by 
an analytically solvable model. This enables us to gain 
physical insight into the essential features of electron 
transport. This model can also serve as the starting 
point for treating one-dimensional electron correlations 
in SACC and their influence on transport. 

The conductance of the system is qualitatively differ- 
ent from that in the case of metal leads. For all electron 
energies corresponding to practically relevant tempera- 
tures and doping levels the junction between the chain 
and the graphene lead is almost perfectly reflecting even 
at strong coupling between the chain and the lead. The 
transmission coefficient of the contact vanishes linearly 
with electron energy as the latter approaches the Fermi 
energy of undoped graphene. This suppression of trans- 
mission results from the vanishing of density of states 
(DoS) in graphene at zero doping. As a result electron 
transport through the interconnect proceeds via narrow 
resonant states in the chain that arise due to strong re- 
flection at the junctions. The width and the position 
of the resonances depend on the length of the intercon- 
nect and the details of its coupling to the leads. The 
shape of low energy resonances is universal but markedly 
different from the Breit-Wigner form. It is dictated by 
the linear energy dependence of the DoS in graphene at 
the point of contact with the chain. This holds even in 
the c ase of leads with zigzag edges, which support edge 
states^ES] Although edge states have a linear disper- 
sion at low energies, their wave functions extend into the 
bulk to distances inversely proportional to the energy and 
give a linear in energy contribution to DoS at the contact 
point. 

Due to the resonant character of transmission the de- 
vice conductance is very sensitive to the number of atoms 
in the chain. In the case of graphene leads the conduc- 
tance difference between chains with odd and even num- 
ber of atoms in the SACC (which was first notecP^ in the 
case of metal leads) becomes much more pronounced and 
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appears as the difference between the on- and off- state 
conductance. 

The high stability of SACC interconnects makes them 
promising building block of atomic scale electronics of 
the future. The resonant character of electron transport 
through them suggests that they can be used as compo- 
nents of atomic scale transistors. Our work is the first 
step towards theoretical understanding of electron trans- 
port in SACC interconnects between graphene leads. 

The paper is organized as follows. In Sec. [II] we qual- 
itatively discuss the essential features of electron trans- 
port in SACC interconnects between graphene leads. In 
Sec. |III| we formulate an analytically solvable model of 
electron dynamics in the system. 
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In Sec 

rive a general formula for the reflection amplitude of the 
junction between the chain and the lead in terms of the 
electron Green function (GF) of the lead. In sec. Ill B 
we evaluate the GF and study the tunneling density of 
states at the junction due to the bulk and edge electron 
states in the lead. In sec. HVI we evaluate the transmission 
coefficient of the device and obtain the universal formula 
for the resonance shape. We discuss our results and ex- 
perimental implications in Sec. |V| 



II. QUALITATIVE DISCUSSION 

Linear molecules with degenerate electron states rep- 
resent an exception to the Jahn- Teller theorem on the 
instability of the symmetric molecular configurations^. 
In this case the degenerate electron states have nonzero 
angular momenta ±A about the molecule axis. Thus the 
matrix element coupling the two states and the corre- 
sponding energy gain arises only in second or higher or- 
der in the vector displacements of the nuclei from the 
symmetry axis. In the linear SACC the two degener- 
ate electron bands have angular momenta ±1 about the 
molecular axis, and their coupling arises only in second 
order in the nuclear displacements. In the Peierls chan- 
nel, which at zero doping corresponds to dimerization, 
the coupling between right and left movers is linear in 
the displacements. Thus the two most likely candidates 
for the SACC structure are the linear chain of equidis- 
tant double bonded carbon atoms, known as cumulcnc 
(•••C = C = C = C- --), and the Peierls-distorted dimcr- 
ized chain, known as polyyne (• • • C— C = C— C • • •). This 
expectation is confirmed by numerical studies. Ab initio 
calculations^ show that because of large quantum fluctu- 
ations of the atomic positions, among all possible spatial 
arrangements of the carbon atoms in the chain only the 
cumulene structure is stable. The outer shell electron 
orbitals in cumulene are sp- hybridized. The two sp hy- 
bridized orbit als form the fully occupied a band. The two 
remaining p y and p z orbitals form two doubly degenerate 
7r bands, which arc half filled making undoped cumulene 
a one-dimensional conductor. Numerical studies^! also 
show that cumulene remains metallic under doping. 

If the SACC attached to the graphene leads is not 



too long the main features of electron transport through 
it can be understood within the single electron picture 
without accounting for electron-electron correlations. 

As the conduction 7r band in graphene leads is formed 
by the p z orbitals, which are perpendicular to the 
graphene plane, only the electrons from the ir z band in 
SACC can propagate into the leads. Thus electron trans- 
port is mediated by a single conducting band in the chain. 

Because of the long mean free path of electrons in 
graphene^, electron motion in the leads may be assumed 
to be ballistic. The device conductance is then deter- 
mined by the elastic electron scattering at the junctions 
between the molecular wire and the leads and backscat- 
tcring in the chain. Backscattering in the SACC due to 
electron-phonon interaction corresponds to emission or 
absorption of phonons with a wavelength of order of the 
interatomic spacing and energy of the order of 10 3 K. As 
a result such processes are exponentially suppresse d 23 * 24 ! 
even at room temperature. Although imperfections in 
the substrate and deviations of the atomic positions from 
the ideal configuration cause some backscattering of elec- 
trons in the wire, we show below that the strongest re- 
flection of the electron wave in the chain occurs at the 
contact with the graphene lead. 

With the aid of the Fermi golden rule the trans- 
mission coefficient of the contact can be estimated as 
T c ~ \ j c \ v w v g- Here j c is the tunneling matrix element 
that couples the last atom in the chain to the lead, v w 
is the local DoS at that atom and v g the DoS an the 
graphene atom that is connected to the chain. The lo- 
cal DoS in the chain is energy independent and can be 
estimated as v w ~ l/7m, where 7^ is the nearest neigh- 
bor hopping integral in the atomic wire. The DoS in 
graphene, on the other hand is strongly energy depen- 
dent and is of the order of v g ~ |e|/7g, where 7 g is the 
nearest neighbor hopping integral in graphene and e is 
the electron energy measured from the Fermi level of un- 
doped graphene^. Thus the transmission coefficient of 
the contact can be estimated as 
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The hopping integrals in graphene and SACC are of the 
same order of magnitude. Therefore at typical doping 
levels, |e| <C 7 ff , the transmission coefficient is small even 
at strong coupling between the chain and the lead, when 
all hopping integrals between nearest neighbor carbon 
atoms are of the same order, j c ~ j w ~ 7 S . 

Neglecting the weak backscattering in the wire we can 
express the energy-dependent transmission coefficient of 
the device, T(e), in terms of the reflection amplitudes 
of the left and right junction between the wire and the 
leads, n/ r (e) = |n/ r (e)| exp(i$j/ r ), as 



T(e) 
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Here (e) | 2 = 1 — |n/ r (e)| 2 are the transmission coef- 
ficients of the junctions, and 4> is the phase accumulated 
by an electron upon returning to the same point in the 
chain after being reflected from both contacts. It can be 
expressed as = 2kJ\f+6i + 6 r , where M is the number of 
atoms in the chain, k is the absolute value of the dimen- 
sionless (measured in units of the inverse lattice spacing 
d of the chain) electron quasimomentum, and Si i r are the 
phases of the reflection amplitudes of the contacts. 

Because at low energies the junctions become strongly 
reflective appreciable transmission through the device in 
this regime occurs only near resonances, where the phase 
cj> equals an integer multiple of 2tt. The energy spacing 
between adjacent resonances is ~ 7c/-A/". 

To obtain a simplified expression for the transmission 
coefficient near a low energy resonance we write the re- 
flection amplitudes of the junctions at low energies as 
|r;/ r (e)| = 1 — c;/ r |e|/7 s , where cu r is a numerical co- 
efficient of the order of unity. This expression follows 
from Eq. ([I]) . Linearizing the energy dependence of the 
phase 4> near the resonance energy eo, 4> — c 4>{ e — e o)/lg, 
where ~ Af is a numerical coefficient, we can write the 
transmission coefficient of the device as, 



T(e) = 



all- 
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Here 7q — (c; + c r ) 2 / (4qc, ) is the transmission coefficient 
at the resonance, and a = c 2 /(4c/c r ) ~ N 2 . 

The width of the resonance is of the order of eo/Af. In 
the case of symmetric contacts, c\ = c r , the device be- 
comes perfectly transmitting on resonance. The shape of 
the resonance is shown in Fig. [4] It is strongly asymmet- 
ric and markedly different from the Brcit-Wigner form 
which arises in the case of metal leads. 

Since low energy transmission through the device pro- 
ceeds via a single resonant state in the chain it is clear 
that Eq. ([3]) holds under very general conditions. The 
assumption that the coupling j c between the chain and 
the lead is energy independent is valid if the resonance 
energy is smaller than the inverse propagation time of 
an electron across the junction. Such resonances always 
exist if the SACC is longer than the junction (for exam- 
ple a small peninsular extending between the graphene 
lead and the chain). As long as the backscattering in 
the chain does not lead to localization the resonant state 
will remain coupled to both leads. The backscattering 
will merely modify its energy and strength of coupling 
to the leads, and can be accounted for by the change of 
parameters in Eq. (3]). Similarly electron-electron inter- 
actions in a finite chain will renormalize the energy and 
the coupling of the resonant state with the leads. 

In the remainder of the paper we present a quantitative 
treatment of simple model of electron transport through 
an cumulene SACC interconnect between graphene leads. 



FIG. 1: Color online. Dependence of the transmission coef- 
ficient T in Eq. ( |30[ ) on e/eo (dashed line) with N = 10 and 
a ~ 1. The solid curve is the Breit-Wigner resonance with 
the same resonance energy eo and width. 



III. SYSTEM AND MODEL 

Consider an ideal cumulene SACC connected to 
graphene leads with perfect zigzag edges, as shown in 

As a first step in the theoretical analysis of the system, 
we assume that the atoms in the SACC are in the ideal 
cumulene configuration. We work in the noninteracting 
electron approximation and describe the electron motion 
in the conducting ir y and tt z bands using the nearest 
neighbor tight binding approximation. More complicated 
band structure of the carbon wire will not significantly 
modify our conclusions. 

Electron transport through the device is fully deter- 
mined by the reflection amplitude of the contact between 
the SACC and the graphene lead. We derive in Sec. |III A| 
a general formula for the reflection amplitude of the con- 
tact in terms of the local DoS in the lead evaluated at 
the atom which is connected to the SACC, Eq. (14). This 



expression is applies for an arbitrary shape of the lead. 
Then in Sec. |IIIB| we specialize to the case of graphene 
leads with zigzag edge. The zigzag edge is likely to be 
formed as a result of an electric failure^ or a tear of a 
graphene strip because it has the least number of dan- 
gling bonds per unit length. We find that the edge states 
present in the case of zigzag edge provide a significant 
contribution to the tunneling DoS at the edge of the 
lead, and thus carry a significant portion of the current 
through the device. 



A. Reflection amplitude of the junction 

Let us first consider a single junction between the 
SACC and a graphene lead. We label the sites in the 
chain by an integer n which enumerates the atoms start- 
ing from the junction, see Fig. [2] The reflection ampli- 
tude of the contact can be found from the retarded Green 
function of the auxiliary system evaluated between two 
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FIG. 2: Schematic picture of the device. The atoms in the 
wire are labeled by n = 1,2,.... The unit vectors of the 
graphene Bravais lattice, ai and a.2, are shown by blue arrows. 
Each unit cell (dashed rhomboid) is labeled by (Ni,N-2). 



points inside the semi-infinite wire using the expression, 

G(n, n') ~ exp[— ik\n — n'\] + r(e) exp[— ik(n + n')]. (4) 

Here k is the absolute value of the energy-dependent 
quasimomentum of the electron in the chain, and the 
Green function is defined in terms of the system Hamil- 
tonian H in the standard way, 

G(n, n') = (n\G\n>) = (n|(e + - iTrV), (5) 

where e + = e + iO. 

We write Hamiltonian of the system as 



H = H w + H g + V, 



(6) 



where H w and H g are respectively the Hamiltonians of 
the semi-infinite wire and the semi-infinite graphene lead, 
and V is the perturbation, which describes electron tun- 
neling between them. 

Introducing the Green function of the unperturbed sys- 
tem , Go = (e+ — H w — Hg)" 1 , we can express the Green 
function of the full system as 



G — Go + GqTGq, 



(7) 



where T is the T-matrix of the junction between the chain 
and the lead given by 



f = (l-VGoyW. 



(8) 



In the nearest neighbor tight binding model the tun- 
neling perturbation V couples only the \n = 1) orbital in 
the chain and a single contact site in the graphene lead, 
which we label as |0). In the 2x2 subspace spanned by 
these states the tunneling perturbation V can be written 
as 



v = ic 



1 

1 o 



(9) 



where 7c is the hopping integral at the contact between 
the chain and the graphene lead. In this case the T- 
matrix depends only on the unperturbed Green function 
within the 2x2 subspace, where it has the form 



Gq 



G 9 (0,0) 







G tt ,(l,l) 



(10) 



Here G g (0, 0) is the Green function of the graphene lead 
at the contact site |0) and G w (l, 1) is the Green function 
of the semi-infinite wire at the site n = 1. 

From Eqs. ^ and Q it is clear that all matrix ele- 
ments of the T-matrix outside the 2x2 subspace vanish. 
Therefore the Green function, Eq. (pjj, within the chain 
can be expressed in terms of the T-matrix of the contact 
as 

G{n,n') =G w (n,n') + G w {n, 1)T(1, l)G w (l, n'). (11) 

Here T(l, 1) = (l|f |1) is the (1, 1) matrix element of the 
T-matrix and 

G w (n,n') = H(e+-i? w )-V) 

is the Green function of an isolated semi-infinite wire. 

We use the tight binding model to describe the electron 
Hamiltonian of the chain, 

oo 

H w = j w ^2 (u w \n)(n\ + \n)(n + l\ + \n + l)(n\), (12) 

n=l 

where j w is the nearest neighbor hopping matrix element 
in the wire, and the on-site energy u w j w describes the 
difference in the work functions between graphene and 
the carbon chain (in our notations the Fermi energy of 
the undoped graphene sheet is set to zero). 

The Green function of the semi-infinite wire can be 
easily determined, see Appendix [A] 

G w (n,n') = — ( e -ik\n'-n\ _ e ~ik(n+n')\ ^gj 



2ij w sin fc V 

Here k is the magnitude of the electron quasimomentum, 
which is related to the energy of the electrons by e = 
u w j w + 2j w cosk. 

With the aid of Eqs. pj, Q and (JlOl we can readily 
express T(l,l) in Eq. (11) in terms of the Green func- 
tion of the graphene lead, G g (0,0). This yields for the 
combined Green function evaluated within the wire, 



G(n,n') = 



ik{r 



2i^f w sin k 



1 
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G 9 (0,0) ^_2ikr, 



1 - 7 e- ifc G 9 (0,0 



Here we introduced a combination of hopping integrals 
in the junction and in the chain, 7 — 7^/7^. 

Comparing the last expression with Eq. Q we obtain 
the reflection amplitude of the junction, 



r(e) = 



l- 7 e' fc G 3 (0,0) 
1 - 7 e" ifc G 9 (0,0)' 



(14) 



Equation ( 14 1 is the main result of this subsection. It 



expresses the reflection amplitude at the contact in terms 
of the Green function of the lead at the contact point with 
the wire, G g (0, 0), and holds for an arbitrary lead. 
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FIG. 3: An infinite graphene plane is separated into two 
halves by adding the perturbation V g , which nullifies the tun- 
neling along the dashed bonds. 



B. Graphene leads with zigzag edges 

We now specialize to the case, in which the graphene 
lead is terminated at the zigzag edge. The zigzag edge 
has the smallest number of broken bonds per unit length. 
It is therefore likely that the gap which appears in the 
graphene strip in the experiments of Ref. [14] is formed 
along this edge. 

We use the nearest neighbor tight binding model to 
describe the electron dynamics in graphene and denote 
the electron 7r-orbitals localized at the atomic sites by 
\A, N) (A sublattice) and \B,N) (B sublattice). Here 
N = (Ni,N 2 ) labels the unit cell with a Bravais lattice 
vector iViai + N 2 a 2 , see Fig. [2] The site \A,N = 0) 
is chosen at the atom which is connected to the carbon 
chain. 



In these notations G g (0,0) in Eq. (14 1 can be ex- 
pressed in terms of the Green function of the semi- infinite 
graphene plane, G g = (e + — H g ), as follows: 



G 9 (0,0) = (AN 



0\G g \A,N 



(15) 



In order to evaluate the Green function of the semi- 
infinite plane G g we start with the infinite graphene plane 
and add the perturbation V g , which nullifies the tunnel- 
ing through the bonds which separate the plane into two 
halves along the zigzag edge, see Fig. [3} 

The Green function of the infinite plane is diagonal in 
the quasimomentum representation due to the transla- 
tion symmetry. We introduce the spinor Bloch functions 
as tf£(N) = exp^JVi + iK 2 N 2 )(ij A (k),^ B (k)), where 
k is the quasimomentum, K± = k • ai and K 2 = k • a 2 
are the projections of the quasimomentum onto ai and 
a 2 , and ipA/B are the wave function amplitudes on the 
A/B sublattices. In the quasimomentum representation 
the (inverse) Green function of an infinite graphene plane 



can be written as matrix in the A/B sublattice space, 



(k) 



"7s 



-iK 2 



1 



-e~ tKl +e lK2 



(16) 

Here 7 g denotes the nearest neighbor hopping integral 
in graphene and we introduced the dimensionless energy 
£+ = t+llg- 

The Dirac points at the corners of the hexagonal Bril- 
louin zone correspond to K\ = K 2 = ±2ir/3. At these 
points the off-diagonal matrix elements in the above 
equation vanish. The Hamiltonian near these points re- 
duces to the familiar Dirac equation with the linear spec- 
trum near the Dirac points as shown in Fig. [4] 

The perturbation V g that cuts the graphene plane into 
two halves is given by (see Fig. |3| , 



V„ = 



lg°N 2 ,N. 



$N 1: 
<5iVi,l 

<5jVi,i<5/v' o 



(17) 



The second matrix in the brackets nullifies electron tun- 
neling between the two halves of the plane, and the first 
matrix describes the on-site potential for the atoms along 
the zigzag edge. This potential is parameterized in our 
model by the dimensionless parameter u g , which is equal 
to the ratio of the on-site potential to the hopping inte- 
gral 7 9 . Because of the diminished number of neighbors 
for the edge atoms the on-site potential is expected to be 
positive and have a magnitude of the order of eV, i.e. of 
the same order as the hopping integral, < u g < 1. 

Due to the symmetry of the problem with respect to 
translations along the edge, (Ni,N 2 ) — > (N 1} N 2 + m) the 
corresponding quasimomentum, K 2 , is conserved. There- 
fore below we use a mixed position/quasimomcntum rep- 
resentation, (Ni,K 2 ). 

In this representation the matrix V g is independent of 
K 2 and has nonzero matrix elements only in the 2x2 
space spanned by the states \ A, N% = 0) and \B, Ni = 1), 
which correspond to the carbon atoms on the opposite 
sides of the divide separating the plane into two halves. 
In this 2x2 subspace V g is given by 



Vg = lg 



-1 



-1 u„ 



(18) 



In this representation the matrix V g is independent of 
K 2 and has nonzero matrix elements only in the 2x2 
space spanned by the states \ A, Ni = 0) and \B, Ni = l), 
which correspond to the carbon atoms on the opposite 
sides of the divide separating the plane into two halves. 
In this 2x2 subspace V g is given by 



Vg = lg 



-1 



-1 u„ 



(19) 



The reflection amplitude of the junction, Eq. (14 1, de 



pends only on the Green function of the semi-infinite 
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graphene inside the same 2x2 subspace. In the mixed 
representation the latter satisfies the equation 

G g {K 2 ) = G , g (K 2 ) + G 0i9 (K 2 )V g G g {K 2 ), (20) 



where the perturbation V g is given by Eq. (19) and 



Go^ g (K 2 ) is the unperturbed Green function inside the 
2x2 subspace (in the mixed representation). The latter 
is evaluated in Appendix [B] and is given by 



la 



_ e _l-C+e 2 - 

'ab\l-C+e 2 +Vab e 



(21) 



where we introduced the notations 



C = 4cos^ 



(l + e) 2 -C, b=(l-e) 2 -C. (22) 




K 2 /n 

FIG. 4: Color online: The solid (black) curve represents the 
intersection of the bulk state spectrum of graphene with the 
K\ — 27r/3 plane that goes through the Dirac point. The 
dashed (red) curve represents the spectrum of the edge states, 
which exist only for 2n/3 < Ki < n. The edge state spectrum 
lies below the bulk state spectrum. 



The branch of V ab in Eq. (21 1 is determined by analytic 



continuation of e from the positive imaginary axis, where 
\/ab takes positive real values. 

Using Eqs. (fl9f, (pOl and d2lj) we obtain 



G g {K 2 ) 



2Vab + (1 - u g )a + (1 + u g )b ( 1 
lg [{l- Ug ya-{l + u g )%] [o 1 



The off-diagonal matrix elements in the above expression 
vanish, as they should due to the absence of tunneling 
between the two half-planes. The (1,1) matrix element 
determines the Green function at the zigzag edge for a 
given quasimomcntum K 2 along the edge. Its imaginary 
part gives the tunneling density of states into the edge for 
a given quasimomentum. It arises from two distinct con- 
tributions of the edge and bulk states, which we discuss 
next. 



1. Tunneling density of states into the zigzag edge 

The tunneling density of states at the zigzag edge of 
graphene is described by the imaginary part of the di- 
agonal matrix elements in the Green function Eq. ( 23 1 . 



Physically, the density of states at the edge contains the 
contributions from the bulk and edge states. The contri- 
bution of the bulk states is described by the imaginary 
part of Vab whereas the contribution of the edge states 
corresponds to the pole at (1 — u g ) 2 a — (1 + u g ) 2 b = 0. 
This condition defines the spectrum of the edge states, 



1 + u 2 g - J(l + u 2 g ) 2 + Au 2 g {2 cos K 2 + l) 



2u a 



(24) 



with cos K 2 < — | . 

This spectrum is plotted in Fig. [4] The inequality 
cos if 2 < — \ reflects the fact that for u g > 0, the edge 
states exist only for e > as shown in the following text. 
For e < 0, the density of states for the edge states vanish, 



or the numerator in Eq. ( 23 ) vanishes together with the 



denominator, eliminating the pole. 

For weak on-site potential at the edge, u g <C 1, the 
edge state spectrum reduces to e = — u g (l + 2cosif2), 
with cos K 2 < — \ ■ In this limit the spectrum and the 
wave functions of the edge states can be understood quite 
easily. In the absence of the on-site potential at the edge, 
u g = 0, these states have wave functions which reside 
only on the A-sublattice and are eigenfunctions of the 
quasimomentum k = (K\,K%). It is easy to see from 



Eq. (24 1 that these states form a degenerate band of zero 



energy states, in agreement with Ref. 1191 From Eq. ( 16 1 



it follows that in order to obey the Schrodinger equation 
in the interior of the lead the quasimomentum of such 
states must satisfy the condition e~ lKl = —1 — e lI<2 . 
Further, since the wave function of these states van- 
ishes on the B-sublattice they remain eigenstates of the 
Hamiltonian even after the plane is separated into two 
halves. The normalizability condition for the edge states 
is Imi^i < 0, implying |1 + e lI<2 \ < 1, which is equiva- 



lent to the inequality below Eq. (24 1. And the amplitude 



of the edge states decay with a factor of 2005(^2/2). 
For weak on-site potential at the edge, u g <§; 1 the edge 
state spectrum may be obtained from the first order 
in perturbation, s{K 2 ) — u g \xJjj\(K 2 )\ 2 , where iPa(K 2 ) 
is the wave function of the edge state at the edge 
atoms. The normalization condition gives 1-0^(^2) | 2 = 
!/E^-oo ex P( 2iV i M i) = -(1 + 2 cos K 2 ). 

The above consideration illustrates that in the pres- 
ence of the on-site potential at the edge the band of edge 
states acquires a finite width of order of the on-site po- 
tential. For strong on-site potential at the edge pertur- 
bation theory is no longer applicable and the spectrum 
of the edge states is given by Eq. (24 1. At zero energy 



the spectrum of these states is linear, which results in the 
finite density of states. It might seem therefore that at 
small energies, e <C 1, the contribution of the edge states 
to the tunneling density of states will be much larger than 
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that of the bulk states. This is not so however because 
at small energies near K 2 — 2tt/3 the wave functions of 
edge states extend into the bulk over many lattice spac- 
ings, so that the local density of such states at the edge 
vanishes linearly with energy As a result, for u g ~ 1 the 
contribution of these states to the tunneling DoS at the 
edge turns out to be of the same order as that of bulk 
states. 

The real space Green function G g (0,0) at the con- 
tact point is obtained by integrating the diagonal ele- 
ment of G g (K 2 )(0,0) in Eq. Q over K 2 : G g (0,0) = 



J Gg{Ki). We write this integral as a contour inte- 
gral over the unit circle of the variable z 2 = exp(iK 2 ). 
Inside the contour the integrand has a simple pole corre- 
sponding to the edge states and a branch corresponding 
to the bulk states. We denote the contribution of the 



pole and the branch cut by G po i e and Gb c respectively 



G s (0,0) = —(G pole + G bc ). 

7g 



A lengthy but straightforward calculation gives, 



(25) 



.|e| (1-«2)0(<O 



% , /3u 2 g ~ 2u g e - e 2 



(26) 



where 0(e) is the step function indicating that the density 
of states due to edge states is present only for e > 0. The 
contribution of the branch cut, Gbc, can be evaluated 
analytically at low energies, 



Gbc 



1 

3w 



iru 



2u g e 



Tog 



3{3u 2 



2u g e — e 2 ) + 3u g — e 



ic- 



(27) 



where c : 



^(l+|-^)w0.22. The imaginary parts 
of both contributions (and with them the tunneling DoS 
at the edge) vanish linearly with energy at small energy. 

2. Reflection coefficient of the junction at low energies 



Substituting the previous Eqs. ( 25 )-( 27 1 into Eq. (14) 



we obtain a simple expression for the reflection coefficient 
of the junction at low energies, |erj <C 1, 



-(e)| 2 = l-77£ 



(28) 



with = -2V4^[(l-^)/^+c]/K + g^-^f ). 
In this regime the electron wave incident from the carbon 
chain into the junction is almost perfectly reflected. 



of order unity defined below Eq. ( 28 I 



Expanding the cosine near a resonance energy e = eq, 
we obtain a simple expression for the transmission coef- 
ficient at small energies, 



1 + a 1 



(30) 



where a = Af 2 n/2 is a dimensionlcss parameter. This 
reproduces the result ^ expected from qualitative con- 
siderations in the case of symmetric coupling. The res- 
onance width is r ~ eo/J\f. The resonance shape is 
strongly asymmetric and markedly different from that 
of the Breit-Wigner resonance, as shown in Fig. [T| 



V. SUMMARY AND DISCUSSION 



IV. DEVICE CONDUCTANCE AT LOW 
ENERGIES: ASYMMETRIC RESONANCES 

The strong reflection at the junction at low energies 
indicates that the transmission coefficient of the whole 
device in Eq. ^ also tends to vanish at small energies 
except in the vicinity of resonances, cos(2A/7c + 2Sq) = 1. 
Substituting |r r /j(e)| = \r(e)\ from Eq. (28) into Eq 
we obtain a simple expression for the transmission coe 
ficicnt of the device at low energies, 



T{e) 



1 



i(2Afk + 2S )} 



(29) 



Here Sq = 2 arctan(7/3w) is the contact scattering phase 
shift at zero energy and k — 2/r/ 2 , where r\ is a number 



We studied electron transport through a single atom 
carbon chain connected to graphene leads. The simplic- 
ity of the hybridization pattern of electron orbitals in 
graphene and carbon chains enabled us to construct an 
analytically solvable model and thereby gain physical in- 
sight into the essential features of electron conduction in 
the system. 

Transmission through the device is dominated by scat- 
tering at the contacts between the chain and the lead. 
For typical temperatures and doping levels in graphene 
the current-carrying electron states have energies much 
smaller than the band width. At these energies the con- 
tact between the chain and the lead becomes almost 
perfectly reflecting. Its reflection amplitude can be ex- 
pressed in terms of the Green's function, G g (0, 0), of the 
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lead at the atomic site connected to the carbon chain, see 



Eq. (14 1. In this equation the parameter 7 describes the 



strength of coupling between the chain and the lead. At 
low electron energies the phase factor e lk may be assumed 
energy independent, as it changes appreciably only at en- 
ergy scales of order of the band width in the wire. In this 
regime the energy dependence of the transmission coeffi- 
cient is dominated by that of the density of states in the 



lead. For graphene leads it becomes linear, see Eqs. (28) 
and ([T}. 

For leads with zigzag edges both the bulk and edge 
states contribute to the DoS at the contact point. Due 
to the difference in the on-site energy between the atoms 
at the edge and in the interior of graphene the band of 
edge states acquires a finite dispersion. The spectrum of 
this band is given by Eq. ( 24 ) and is plotted in Fig. [3] Al- 
though the edge state spectrum is linear at small energies 
its contribution to the local DoS at the edge is not con- 
stant, but rather is linear in the electron energy, ~ e#(e)- 
This occurs because the edge state wave functions extend 
into the bulk to distances which are inversely propor- 



tional e, as explained in Sec. IIIB 1 As the difference in 



the on-site potential between the atoms at the edge and 
in the interior of the lead is of the same order as the band 
width the contribution of edge states to the DoS is of the 
same order as that of the bulk states. Therefore a sub- 
stantial part of the current through the carbon chain is 
propagated into the lead by the edge states. The energy 
dependence of the reflection coefficient of the junction is 
described by Eq. (28 1. 



The interference between reflection amplitudes of the 
left and right junctions gives rise to the transmission co- 
efficient of the device described by Eq. ( 29 1 . Due to the 



nearly perfect reflection at the contact the energy depen- 
dence of the transmission coefficient of the interconnect 
has resonant character. Near the resonance the trans- 
mission coefficient is described by a simple expression, 
Eq. pi 



Our main conclusions, namely the linear energy de- 
pendence of the transmission coefficient of the junction 
between the chain and the lead and the shape of the 



resonance in Eq. ( 30 1 do not depend on many of the sim- 



plifying assumptions of our model. 

The linear energy dependence of the junction transmis- 
sion coefficient holds if the coupling between the chain 
and the lead is energy independent. This assumption is 
valid as long as the electron energy is smaller than the 
inverse propagation time across the contact and holds for 
more complicated junctions, e.g. a small peninsular con- 
necting the chain to the lead. In this case Eqs. (14) and 



(29 1 will still hold, provided 7 and r\ are replaced by the 



appropriate parameters describing the coupling strength 
between the chain and the lead at low energies. Similarly, 
Eq. ( 30 ) will also hold provided the resonance energy and 



the parameter a are chosen appropriately. The general- 
ization of the resonance shape to the case of asymmetric 
contacts is given by Eq. ([3|. 

The resonant character of transmission will be pre- 



served even in the presence of the Coulomb interaction in 
the wire, as long as wire is short enough so that the one- 
dimensional correlation effects can be neglected. Such a 
wire will act as a molecule with a single resonant level 
participating in transport. For longer wires the one- 
dimensional correlations need to be taken into account. 
In this respect the Umklapp processes and the formation 
of Friedel oscillations near the contact points are espe- 
cially important. The study of these effects is left for 
future work. 
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APPENDIX A: GREEN FUNCTION OF A 
SEMI-INFINITE WIRE 

We construct the Green function of the semi-infinite 
wire from that of the infinite wire by adding a perturba- 
tion that nullifies the hopping between the two halves. 

The retarded Green function of the infinite wire in k 
space is diagonal and given by 



with 



G , w (k) = (e-e k + i0)- 



£fc = u w -f w + 2j w cos k. 



(Al) 



(A2) 



The real space Green function is obtained by integrating 
over k as 

G (n,n') = i- I dkG . w (k)e^ n - n "> (A3) 



which gives 

G . w (n,ri) = 



2ij w sin k 



exp(— ik\n — n |) 



(A4) 



where k is the magnitude of the electron quasimomcntum 
related to the energy by e = u w ~f w + 2j w cos k. 

The perturbation V w which cuts the wire to two halves 
has non- vanishing matrix elements only in the 2x2 sub- 
space spanned by the orbitals with n = and n = 1, 
where it is given by 



V w (0,0) V W (0,1) 
V W (1,0) V W {1,1) 



7™ 

lw 



(A5) 



The T matrix defined in Eq. ^ is also nonvanishing only 
in the 2x2 subspace and can be expressed solely in terms 
of the matrix elements of Gq. w in the 2x2 space, 



Go,40,0) G ,»(0,1) 
Go,mi(1,0) Go,m,(1,1) J 2i^ w sink \ •• 



1 e 

— ik ^ 



■/ k 



(A6) 

Using Eqs. (|A4]), flAS) , ( |A6| , Q and (|8) we obtain the 
Green function of the semi-infinite wire, Eq. (13). 
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APPENDIX B: DERIVATION OF EQ. (21) 



calculation gives 



In this appendix we derive the expression for the un- 
perturbed graphcne Green function within the 2x2 sub- 
space spanned by the rows of atoms on the opposite sides 
of the dashed links in Fig. [3] 

The unperturbed Green function in the quasimomen- 
tum representation is obtained by inverting the matrix 
in Eq. (fl6|, 



Go, s (k) = - 



— e 

-iJCi 



-iK x 



l + e 



-iK 2 



where D — j g [e 2 — (1- 

In the mixed representation the Green function 
Go t g{K<i) in the 2x2 subspace of states \A,N\ = 0) 
and \B,Ni — 1), can be obtained by the inverse Fourier 
transform of Go, 9 (k) with respect to K±. An elementary 



and 



(A,0\G JK 2 )\A,0) = (B,l\G , g (K 2 )\B,l) 
dK\ e e 



2tt D 



Is 



ab 



(Bl) 



(A,0\G , g (K 2 )\B,l) = (B,l\G , g (K 2 )\A,oy 



dKi l+e 



2tt 



IgVab 



D 

l-C + e 2 



'ab 



(B2) 



where a, b and C are defined in Eq. (22 1 



Combining the above matrix elements into one 2x2 
matrix we arrive at Eq. (211. 
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It might seem that the presence of edge states might give 
rise to an energy-independent contribution to the tunnel- 
ing DoS at the point. This is not so however. The wave 
functions of edge states at low energies extend into the 
bulk to distances inversely proportional to the electron en- 
ergy. As a result their contribution to the local DoS at the 
contact is also linear in e. It is analyzed in detail in section 
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